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Abstract 



We provide compeUing evidence that a previously introduced model of non-perturba- 
tive 2d Lorentzian quantum gravity exhibits (two-dimensional) flat-space behaviour 
when coupled to Ising spins. The evidence comes from both a high-temperature 
expansion and from Monte Carlo simulations of the combined gravity-matter sys- 
tem. This weak-coupling behaviour lends further support to the conclusion that 
the Lorentzian model is a genuine alternative to Liouville quantum gravity in two 
dimensions, with a different, and much 'smoother' critical behaviour. 
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1 Introduction 



At the end of the twentieth century, the non-perturbative quantization of gravity 
remains an elusive goal for theoretical researchers. There is not even a consensus 
on how the problem should best be tackled. For example, considering pure-gravity 
approaches, we have on the one hand Euclidean path-integral methods, which are 
close to usual formulations of (non-generally covariant) quantum field theories and 
well-suited for numerical simulations. On the other hand, in canonical quantization 
approaches it is - at least in principle - easier to address questions about the be- 
haviour of spatial three-geometries, but the complicated structure of the constraints 
tends to lead to computational difficulties.^ Unfortunately, very little is known 
about the relation between the covariant and canonical approaches. In part this is 
due to the 'signature problem' of the path-integral formulations: the sum over all 
space-time geometries is usually taken over Riemannian, and not over the physical 
Lorentzian (pseudo-Riemannian) four-metrics modulo diffeomorphisms. The prob- 
lem of how to relate the two sectors by an appropriately generalized Wick rotation 
remains unresolved. 

Our aim is to investigate the possible consequences of taking the Lorentzian 
structure seriously within a path-integral approach. In order to gauge the difficulties 
this involves and to circumvent technical problems, we first addressed the issue in 
two space-time dimensions, where there already exists a well-understood theory of 
(Euclidean) quantum gravity, namely, Liouville gravity. In we proposed a new, 
Lorentzian model of 2d quantum gravity, obtained by taking the continuum limit of 
a state sum of dynamically triangulated two-geometries. The Lorentzian aspects of 
the model were two-fold: firstly, the sum was taken only over those two-geometries 
which are generated by evolving a one-dimensional spatial slice and which allow 
for the introduction of a causal structure. Secondly, the Lorentzian propagator 
was obtained by a suitable analytic continuation in the coupling constant. The 
first aspect turned out to be the crucial one, leading to a continuum theory of 2d 
quantum gravity inequivalent to the usual Liouville gravity. This was shown in 

where both the loop-to-loop propagator and various geometric properties of the 
model were calculated explicitly. The Hausdorff dimension of the quantum geometry 
is (ii^ = 2, and points to a much smoother behaviour than that of the Euclidean case 
(where (i// = 4). 

However, we must emphasize that (1^ = 2 does not imply a flat geometry. The 
model of Lorentzian gravity defined in allows for arbitrarily large fiuctuations 
of the spatial volume from one time-slice to the next. This is illustrated by Fig. 
ID, which shows a typical surface generated by the Monte Carlo simulations, to be 
described in Sec. |[ The length of the compact spatial slice fiuctuates strongly with 
time (pointing along the vertical axis). Using the results of [|I|, one easily derives 

^Alternatively, one could embed quantum gravity in a larger, unified theory like string theory 
or (the as yet non-existent) M-theory. However, these are still far from giving us any detailed 
information about the quantum gravity sector. 
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Figure 1: A typical discrete history of pure Lorentzian gravity with volume N=1024. 



that in the thermodynamic limit and for large times the average spatial volume L 
and fluctuations around L behave like 



respectively, for a given cosmological constant A. This demonstrates that even in 
the continuum limit the fluctuations are large, and of the same order of magnitude 
as the spatial volume itself. 

We managed in [jl| to further trace the difference between the two quantum the- 
ories to the presence or absence of so-called baby universes. These are outgrowths 
of the two-geometry giving it the structure of branchings-over-branchings, which are 
known to dominate the typical geometry contributing to the Euclidean state sum. 
On the other hand, in the Lorentzian state sum, one can suppress the formation 
of such branchings with respect to the preferred spatial slicing (which is not present 
in the Euclidean picture, where no directions are distinguished). There is also a 
physical motivation for suppressing the generation of baby universes, since the asso- 
ciated (discrete) geometries can usually not be embedded isometrically in a smooth 
Lorentzian space-time. If nevertheless one did decide to generalize the evolution 
rules of the Lorentzian model to allow for such branchings (and keep only a weaker 
notion of causality, c.f. |ll]), one would rederive the usual Euclidean Liouville re- 
sults. In what follows, when talking about 'the Lorentzian model', we will mean the 
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unmodified model without branching baby universes, i.e. the model of 2d quantum 
gravity that does not lie in the same universality class as Liouville gravity. We also 
would like to point out that from the point of view of a canonical quantization the 
Lorentzian model is much more natural. The inclusion of topology changes of space 
into a canonical scheme would require a so-called third quantization of geometry. 

In Liouville gravity, matter couples strongly to geometry, perhaps even too 
strongly in the sense that the combined system becomes inconsistent when the cen- 
tral charge of the (conformal) matter exceeds one. Arguments have been presented 
which link the strong deformation of geometry to the creation of baby universes 
0. It is therefore conceivable that the Lorentzian model of gravity - where baby 
universes are absent - has a weaker and less pathological coupling to matter. 

In order to understand the behaviour of the combined gravity-matter system, 
we are considering here the coupling of the gravitational model of reference ^ to 
an Ising model of spin-i, with nearest-neighbour interaction J2<ij><^i<^j between 
its spins cTj = ±1. In [Q] we made a careful analysis of the implications of the 
Lorentzian signature for the sum over space-time metrics. The most straightforward 
way of obtaining the continuum limit consisted in performing the calculations in 
the discretized model with purely imaginary coupling (corresponding to Euclidean 
signature) and only afterwards 'rotating back' to the Lorentzian sector. Moreover, 
it turned out that certain simple properties, like the fractal dimension of space-time, 
were independent of the analytic continuation. 

We will apply the same philosophy in the present context by analyzing the Ising 
model coupled to 2d gravity with coupling constants corresponding to the Euclidean 
signature sector. Nevertheless we will continue to talk about 'Lorentzian' gravity 
coupled to matter, because the choice of two-dimensional Euclidean geometries con- 
tributing to the path integral is dictated by the requirement that after the rotation 
to Lorentzian signature they should be causal and non-singular. The 'matter observ- 
ables' we will consider are the critical exponents for the Ising model, characterizing 
the underlying c= 1/2 fermionic continuum model coupled to gravity, which are 
not expected to change under the rotation to Lorentzian signature. In order to 
determine the universality class of the interaction between matter and gravity it is 
therefore convenient to work entirely within the Euclidean sector of our Lorentzian 
gravity model. 

For fixed regular two-dimensional lattices, and in the absence of an external 
magnetic field, the Ising model can be solved exactly in a variety of ways (see, 
for example, ^, ^). The partition function (for the square lattice) was found 
by Onsager. Its critical behaviour is characterized by a logarithmic singularity of 
the specific heat and the critical exponents near the Curie temperature, a = 0, 
(3 = 0.125, and 7 = 1.75, for the specific heat, the spontaneous magnetization and 
the susceptibility respectively. 

For the case of the usual Euclidean 2d gravity, described by an ensemble of planar 
random surfaces, coupling to Ising spins was first considered in 0, where an exact 
solution was obtained with the help of matrix model methods. It could be shown 
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that in the presence of gravity, the matter behaviour is 'softened' to a third-order 
phase transition, characterized by critical exponents a = —l, (3 = 0.5, and 7 = 2 [0. 
On the other hand, the geometry is 'roughened', as exemphfied by the increase from 
— 1/2 (pure 2d Liouville gravity) to —1/3 of the entropy exponent ■y string for baby 
universes on manifolds of spherical topology. 

It is not entirely straightforward to apply the methods used to obtain these ex- 
act solutions to the Lorentzian gravity model. For example, one can write down an 
expression for the transfer matrix generalizing that of the Onsager solution by im- 
posing a length cutoff Iq on the length of spatial slices. However, a major stumbling 
block to understanding the behaviour of its eigenvalues as /q 00 is the fact that 
as a consequence of the gravitational degrees of freedom, transitions between spa- 
tial slices of different length are allowed. This makes the use of Fourier transforms 
problematic, which are an essential ingredient of this and other algebraic solution 
schemes. Moreover, the Hilbert space dimension for the discrete, finite model is 
given by E1Li2', which grows rapidly with Iq. 

In the absence of an analytic exact solution^], one way to try to extract infor- 
mation about the matter-coupled model is by performing a series expansion of the 
partition function Z at high or low temperature, or of suitable derivatives of Z. For 
flat, regular lattice geometries, these have been studied extensively since the early 
days of the Ising model. It is well-known that the high-T expansion, in particular, 
that of the magnetic susceptibility x ^^ro field is well-suited for obtaining infor- 
mation about the critical behaviour of the theory. We will show that the same is 
true for the coupled gravity-Ising model, after taking into account some peculiarities 
to do with the fact that we have an ensemble of fluctuating geometries instead of a 
fixed lattice. In the limit of large lattice size A^, there is a well-defined expansion in 
terms of m := tanh/3, where the coupling (3 is proportional to the inverse temper- 
ature, whose coefficients can be determined by diagrammatic techniques. Given a 
plausible ansatz for the singularity structure of the thermodynamic functions, one 
can then extract estimates for the critical point and critical exponents from the first 
few terms of such an expansion. These results are corroborated by performing a 
Monte Carlo simulation of Lorentzian gravity coupled to the Ising model. Apart 
from being in good agreement with the high-T expansion, the simulations also allow 
us to measure the quantum geometrical properties of the model. 



matrix model of Lorentzian gravity coupled to Ising spin has been formulated recently. Its 
analysis is the subject of a forthcoming publication pO[ |. 
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2 The high-T expansion 



Recall the usual high-T expansion of the Ising model on a fixed lattice of volume 
A^, with partition function 

Z(/?,iV)= e " , /3 = 77^, (2) 

where the sum is taken over all possible spin configurations, and J > denotes the 
ferromagnetic Ising coupling. We will only consider the case of vanishing magnetic 
field, H = 0. The Ising spins are located at the lattice vertices, labelled by i,j G 
1 . . .V. A convenient expansion parameter at high temperature is m := tanh f3, which 
we can use to re-express 

e^"'"^ = {l+u a.Gj) cosh (3) 
Substituting (|) into (||), the partition function becomes 

N) = (cosh/3)^ E [l + « E ^i^J + E E (^.^.OK^O + • • •] (4) 

{crj <ij> <ij> <kl> 

=: 2'^(cosh/3)^(l + 5]^]„M"), (5) 

n>l 

with V denoting the number of vertices and s the number of nearest-neighbour pairs 
(i.e. the number of lattice links). Note that the terms ~ in eq. (|^) are only non- 
vanishing if every ai in ai^cxi^ • • • cTj^ appears an even number of times. Representing 
spin pairs (o"j(Tj) by drawing a link between cxj and aj on the lattice, this is equivalent 
to the following statement: no n- vanishing contributions to Qn in eq. (|^) correspond 
to figures of lattice links which are closed polygons, with an even number of links 
meeting at each vertex. The coefficient Qn simply counts the number of such figures 
at order n that can be put down on a given lattice, and will depend on the lattice 
geometry (triangular, square, etc.). It is a polynomial in the variable A^. 

Because of the extensive nature of the free energy F{N) = —kT In Z{N) , we must 
have that (1 + ^nu"') ~ e^*^'"-' in the thermodynamic limit — oo, and we can 
therefore write for the partition function per unit volume 

\nZ{(3) := llnZ(/?, AT) = ^ cosh/3+ ^ ln2 + (6) 

where cu^^ is obtained by taking the term linear in A^ in f2„ and setting A^ = 1. 
Note that both connected and disconnected graphs contribute to u!^\ A simi- 
lar relation can be obtained for the magnetic susceptibility at zero field, x{^) = 
kT-^j2 In Z{N)\h=o- At high temperature, the susceptibility per unit volume can be 
expressed as 

X = kT{l + Y.^f^u-). (7) 

n>l 



6 



The coefficients tu^^-* are the exact analogues of cj^*^-* in eq. (j^), where the counting 
now refers to polygon graphs with two odd vertices (vertices with an odd number 
of incoming links), and all other vertices even (c.f. 0, but beware of the difference 
in notation for the number of vertices). 

Since we are primarily interested in the bulk behaviour of the gravity-matter 
system, we will in the following for simplicity choose the boundary conditions to be 
periodic. That is, we will identify the top and bottom spatial slices of the cylindrical 
histories introduced in 0]. Clearly this is not going to affect the local properties 
of the model. As above, we will denote the discrete volume, i.e. the number of 
triangles of a given two-dimensional geometry (with torus topology), by N. It 
follows immediately that such a geometry contains N time-like links, N/2 space-like 
links, N/2 vertices and 3N/2 nearest-neighbour pairs. 

In quantum gravity the volume becomes a dynamical variable. For fixed 
topology, the only coupling constant appearing in the action of pure 2d quantum 
gravity is the cosmological constant, multiplying the volume term. The partition 
function of the Ising model coupled to 2d Lorentzian quantum gravity is given by 

GiX,t,P) = Y: e-'^'-ZTiP) = E e-'^" E e^-?--'"'"', (8) 

where the sum is taken over all triangulations T with the topology of a torus and 
t time-slices, Nt is the number of triangles in T, and Zt{(3) the Ising partition 
function (|^) defined on T. Fortunately, the summation over volumes in eq. (§) does 
not lead to additional complications in the analysis of the thermodynamic properties 
of the spin system, since the state sums for fixed and fluctuating volume are simply 
related by a Laplace transformation. Rewrite relation (§) as 

G(A,t,/?) = Ee-'^^(AA^,t):=Ee"'^ E ^t(/?), (9) 

where T/v.t denotes the toroidal triangulations of volume and length t in the time 
direction. Analogous to eq. (1), we expect the matter part /(/?) of the free energy 
density in the gravitational ensemble to behave in the thermodynamic limit (A^ oo 
and t oc VN) likeQ 

Z(/?,Ar,t)^e("-^^('3))^+''(^). (10) 

(For simplicity, we have set the ferromagnetic coupling to J = 1.) We can now 
reexpress eq. (|) as 

G(A, t, /3) = E e(^=(^)-^)^+°(^) , A,(/3) = A, - /3/(/3), (11) 

N 

where Ac = XdP = 0) = ln2 is the critical cosmological constant of pure gravity, 
which was determined in |jl|. Interesting limiting cases are — > where —I3f{(3) = 

^Note that with the conventions used in definition (^), the ground state energy is —3PN/2 and 
the free energy density f{(3) is negative. 
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I In 2, reflecting the factor 2^ in eq. (|) (eacli spin lias two states), and the strong 
coupling region /? — > oo where —Pf{P) 3/3/2 (only the ground state of all spins 
aligned contributes to the state sum (^). The term proportional to the pure gravity 
cosmological constant Ac appearing together with the free energy in ([T0|) has its origin 
in the sum over all triangulations, 

A calculation of Z{(3, N, t) not only determines the thermodynamic properties of the 
spin system in the presence of gravity, but at the same time describes gravitational 
aspects of the coupled system, for example, the critical cosmological constant XdP). 
Conversely, knowledge of Ac(/3) determines the spin partition function in the infinite 
volume limit. - The analogue of the high-T expansion in the presence of gravity 
is given by 

Z(AiV,t) = (cosh/3)^ 2T ^ (l + ^17„(T)n"). (13) 

TeTjv.t n>l 

We may reexpress the critical cosmological constant of the combined system as 

Ac(/?) = A, + ^ In cosh /3 + 1 In 2 + f{u) , (14) 
where f{u) is defined in the thermodjTiamic limit by 

E (1 + E ^n{T)u^) 

^^^^'^ ""^^ = e^fi'^). (15) 

Terjv,t 

The coefficients Q of the power series now depend on the triangulation T. When 
counting diagrams of a given type and order n, we must keep in mind that the 
vertex neighbourhoods do not look all the same, as they do in the case of a regular 
lattice, but that the distribution of coordination numbers (numbers of links meeting 
at a vertex) is subject to a probability distribution. The coefficients in the high-T 
expansion therefore count the average occurrence of a certain diagram type in the 
ensemble of triangulations of a fixed volume N, for large A^. 

Starting to evaluate the series in (|T3D order by order, one immediately notices a 
qualitative difference from the regular case. If we had considered a regular triangular 
lattice (coordination number 6), the first non-trivial contribution to the counting 
of even diagrams would have appeared at tt, = 3, where one obtains Q^lN) = N, 
coming from closed triangle graphs. However, when looking at all two-dimensional 
random lattices contributing to the sum over geometries in the gravity case, there 
are geometries which have one or several 'pinches'. A pinch is a spatial slice of 
minimal length 1 = 1, which consists of a single link and a single vertex (see Fig. 0). 
Pinches occur even if the total volume of the two-geometry is kept fixed, since in the 
presence of gravity the length of spatial slices is a fiuctuating dynamical variable. 
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Figure 2: A two-dimensional geometry with a 'pinch' of length 1. 

For the gravitationally coupled Ising model, the lowest-order contribution to 
the power series in u in ([l3[) occurs therefore already at order n = 1. Clearly, 
such pinching contributions will be present at all orders, in both connected and 
disconnected diagrams, on top of the ordinary 'bulk' contributions, coming from 
diagrams which do not wind around the spatial direction of the torus in a non- 
trivial way. The former have no analogue on regular lattices. 

Fortunately, it turns out that the pinch contributions are irrelevant, in the sense 
that they contribute at a lower order of A^, whereas the bulk contributions in f2„ go 
like N'', k > 1. This can be seen most easily by considering the Laplace-transformed 
partition function. Let us begin by evaluating the zeroth-order term of eq. (|13|), 

G(A,t) = Ee-'^ E l:=Ee-''^(cosh/5)^2f h (16) 

where for notational brevity we have defined an 'effective cosmological constant' 
A = A — I In cosh/? — |ln2, in accordance with eq. The left-hand side of (|16|) 

can be computed as 

given the explicit form of the propagator derived in ref. to which we also refer 
for details of notation. The term proportional to in the Laplace transform of (p!3D 
is 

r^u^: Ee"^"^ E l = G'(A,t)fir"(A), (18) 

where the second summation is over triangulations with a single 'pinch' of spatial 
length / = 1. To arrive at the last expression on the right-hand side, the factor 
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G(A,t) has been pulled out. In terms of quantities derived in the normalized 
coefficient is most easily computed as 

EG'^{x,l = l;i)G'^{l = l,y;t-i) 

= — 777 T^ • (19) 

G-^{x,y;t) x=y=o 

We are interested in the behaviour of this expression in the thermodynamic limit, 
which is tantamount to letting the cosmological constant approach its critical value, 
A — s> Ac- In this limit, ([T9|) yields simply a constant, 2. This is a general 



feature of configurations with one or several pinches. For example, generalizing to 
geometries with a single pinch of length / gives a coefficient 21 in the large-volume 
limit. As an example of a more complicated configuration, the normalized coefficient 
for histories with one pinch of length li and a second one of length I2 becomes in 
this limitQ 



~ u 



min(ii,/2) — 1 



3 E \ \L u'.ur (20) 



k=0 



{h - k - i)\{i2 - k - ly.ki' 



By contrast, let us now calculate the first bulk contribution, which occurs at 
order u^. The contribution to ^3 is simply A^, from counting the number of triangle 
graphs in the 2d geometry. Taking the Laplace transform, we obtain 

E^-'^'N ^ -^G{X,t) ^ {N) G(A,t). (21) 

Evaluating the expectation value of A^ in the continuum limit, one finds 

(AT) = -G{X,t)-' ^GCX,t) (22) 

aA 



4 (1 _ e-^TVA _ T^(l - e-2^v^)) t large 4T 



a2A(l - e-2^v^) a^VX' 



(23) 



(We are using the notation of [|T|, with T and A the continuum length of the two- 
geometry in 'time'-direction and the renormalized cosmological constant.) This di- 
verges exactly the way one would expect from a volume term. It reiterates the con- 
clusion of [0] that all macroscopic metric variables scale canonically in the Lorentzian 
gravity model. 

^Let us take the opportunity to correct some misprints in equation (29) of which has been 
used in deriving formula (EG). The correct equation reads 



l2BlA\'+'- ^ {h-k-l)\(h-k-l)\k\\ Bf 



where F, At , Bt and Ct are defined in |^ 
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i-1 i i+1 



Figure 3: The triangles contributing to the weight at the vertex i. 

We conclude that in the thermodynamic limit, pinching terms will be suppressed 
since their number is proportional to N^, whereas the (connected) bulk diagrams 
behave like ~ A^^. For large A^, the pinch contributions must therefore factorize 
according to 

(1 + ^ f2„(T)M") = {1 + N^Y. PmiT)un (l + NYl ^n{T)u^ + 0{N^)). (24) 

n>l m>l n>l 

Taking the logarithm, we see that the sum (1 + J2Pm{T)u"^) ~ will only con- 
tribute a constant term to the free energy, which does not affect the universal be- 
haviour of the model. We will make no attempt to calculate it explicitly. Similar 
considerations apply to the high-T expansion of the magnetic susceptibility in the 
presence of gravity. The pinch contributions factorize, and we will only need to com- 
pute the multiplicity lvI^^ of bulk polygon graphs with two odd vertices per triangle 
in 

X~(1 + E4M- (25) 

n>l 

Our next step will be to derive the probability distribution of the coordination 
numbers in the Lorentzian gravity model, in the thermodynamic limit as the cos- 
mological constant A Ac = In 2. Recall that when generating an interpolating 
space-time between an initial and a final spatial geometry, the geometry of each 
space-time 'sandwich' with At = 1 is independent of the previous one in the sense 
that there are no local constraints on how the numbers fc, > 1 of time-like future- 
pointing links can be chosen at each vertex i Having reached a spatial slice at 
time t, we can generate the space-time between t and t + 1 proceeding from 'left to 
right'. To each vertex i at time t we associate ki time-like links (ending at vertices of 
the subsequent spatial slice at t + 1) and the space-like link to the right of the vertex. 
There are therefore exactly ki triangles associated with the vertex i, contributing 
with a weight factor e~''^ to the action, as illustrated in Fig. |^. Since the assignment 
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of the order ki of outgoing time-like links to the vertex is completely independent 
of the ^-assignments of other vertices, the probability distribution for k outgoing 
future-directed links is given by 

pxik) = e''\e^ - 1). (26) 

Strictly speaking, the argument leading to eq. (pB]) is only correct in the continuum 
limit in which 'extreme pinching' to vanishing spatial length / = does not occur 
(for off-critical A, relation ( pB]) must be modified to account for the fact that moves 
changing the torus topology are forbidden). Fortunately this is the only case we are 
interested in, and the final probability distribution is therefore obtained by setting 
e"'*' = 1/2 in (^), yielding 

p{k)^p,M = Yk- (27) 

For reasons of symmetry, the distribution of incoming time-like links at i (originating 
at the slice at t—1) is of course identical. Given relation (^), we can now compute 
the probability distribution p{j) of the vertex order, i.e. of the total number j of 
links meeting at a vertex (incoming and outgoing time-like and space-like links), 

PU) = J> 4. (28) 

With the distribution ( pTf ) in hand, we can now embark on the actual counting 
of diagrams contributing to the susceptibility coefficients tu^^-' in (P5|). We will only 
quote the results up to order n = 5. Further details of the counting procedure will 
appear elsewhere. The average numbers of diagrams per triangle (i.e. per unit 
volume) are listed in the table below. 



n 


open 


closed 


disconnected 


total 


1 


3 
2 








3 
2 


2 


si 








8i 


3 


43i 








43i 


4 


214| 


14 


-17 


211| 


5 


1038i 


134-^ 


-174— 

' ^18 


997i 



Open graphs are connected graphs without any self-intersections. Closed graphs are 
connected graphs which are not open. The disconnected graphs consist of two or 
more components and contribute with a minus sign. 

In order to double-check our results at order 4 and 5, where the counting becomes 
slightly involved, we have performed a numerical check on the coefficients tu^^^ listed 
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above. This was done by computer-generating histories of length At ~ 100, with 
an initial spatial slice of length Al = 200, and counting diagrams of a given type. 
The results are given in the table below and in very good agreement with the exact 
calculation. They are based on a total of ~ 3 x 10^ vertices at order 4 and ~ 
9 X 10^ vertices at order 5. We have not listed the counting of disconnected graphs 
separately, since it follows closely the counting of closed connected graphs. 



n 


open 


closed 


4 


214.642 ±0.179 


13.996 ±0.007 


5 


1037.770 ±0.751 


134.197 ±0.098 



2.1 Evaluation of results 

In order to evaluate the results from the high-T expansion, we assume a simple 
behaviour of the susceptibility of the form 

x{u)r^A(l--y' + B (29) 



Uc 

near the critical point u^, with analytic functions A and B. Using the ratio method 
(see, for example, ||lll), we have fitted the susceptibility coefficients to 

.,. = i|^ = l(n-lzl). (30) 

Plotting the ratios r„ linearly against 1/n for n G 1 . . . nmax, we have extracted the 
following estimates for the critical point Uc and the critical susceptibility exponent 
T 



^max 


critical point 


critical exponent 


3 


= 0.2488 


7 = 1.820 


4 


Uc = 0.2462 


7 = 1.789 


5 


Uc = 0.2458 


7 = 1.783 



The estimates for the critical exponent should be compared to the exact values 
for 7 for the Ising model on a fixed, regular lattice and on dynamically triangulated 
lattices (Ising spins coupled to Euclidean quantum gravity), which are = 1-75 



13 



and 7*^* = 2 respectively. The data from the high-T expansion clearly favour 7 = 1.75 
in our model. Indeed, the estimates for 7 are remarkably close to this value, given 
that we are working only up to order 5 in the expansion parameter u = tanh/3. 
The conclusion that the critical exponents of the Ising model coupled to Lorentzian 
quantum gravity coincide with those found on regular lattices is also supported by 
the Monte Carlo simulations we have performed. 

However, before turning to a detailed description of the simulations we would like 
to illustrate how well the high-T expansion works even at this rather low order. We 
will compare the /3-dependent cosmological constant XdP) defined in eq. (^^, which 



can be measured directly in the Monte Carlo simulation, with the same quantity 
obtained from the high-T expansion. Recall that in the thermodynamic limit Ac(/9) 
is essentially given by the spin free energy, eq. (|llD, which can be computed in the 
small-/? expansion. We have determined the density f{u), defined in eq. (|T5|) , by 
counting closed polygon graphs in the high-T expansion up to order 6. Inserting 
this into formula f|l3) leads to 



A^'s^-^(/3) = A, + i In 2 + ^ In cosh/3 + + ^ + fu' + '^u'. (31) 

In Fig. ^ we show the data points for Ac(/5) — 3/?/2 as measured by the Monte 
Carlo simulation]^. Since Ac = In 2 in pure gravity, the data should approach | In 2 
for /? — > and Ac = In 2 for f3 — >■ 00, both of which are well satisfied. In order 
to quantify the effect of the w-expansion, we have plotted both the zeroth-order 
expression Fi{P) = Ac + |ln2 + | In cosh/? — |/?, and the improved sixth-order 
expression F2{P) = Aj}^^*^~"'"(/5) — |/3. The latter agrees well with the measured 
Monte Carlo values right up to the neighbourhood of the critical Ising coupling (3^- 
At the critical point jSc the measured function Ac(/?) exhibits a cusp. This reflects 
the singular part contained in XdP) which of course cannot be captured by simply 
plotting the analytic function ([T^). 



3 The Monte Carlo simulation 

Monte Carlo simulations have been used successfully in the study of Euclidean 2d 
quantum gravity. The formalism known as 'dynamical triangulations' provides a 
regularization of the functional integral well-suited for such simulations, allowing in 
addition for a straightforward matter coupling of Gaussian fields as well as of spin 
degrees of freedom. Extensive computer simulations of the combined gravity-matter 
systems have been performed, leading to results in perfect agreement with exact 
results derived from Liouville theory and matrix model calculations. 

^The subtraction of 3/3/2 has been performed to ensure a finite limit as /3 — > 00. It corresponds 
to using the action f3J2<ij>i^'i''^j ~ 1) ™ '^'i- (H)' whose ground state has energy zero rather than 
-3/37V/2. 
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Figure 4: The critical cosmological constant as a function of the Ising couphng (3, 
as measured by Monte Carlo simulations {t — 32, N — 20AS), and compared to the 
corresponding high-T expansions -Fi(/3) and F2{/3) at order and order 6. 
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Figure 5: The move used in the Monte Carlo updating of the geometry. 



The Lorentzian model resembles the dynamically triangulated model in that its 
dynamics is associated with the fluctuating connectivity of the triangulations con- 
tributing to the path integral. This allows us to take over many of the techniques 
from the computer simulations of the dynamically triangulated models. We must 
specify the update of both the geometry and the matter fields, the latter being stan- 
dard: for a given triangulation we update the spin configurations by the same spin 
cluster algorithms used for dynamical triangulations. This presents no problems 
since our configurations form a subset of the full set of dynamical triangulations 
used in Euclidean quantum gravity (on the torus). During the update of geometry, 
we want to keep the number of time-slices fixed while allowing any space-like fluc- 
tuations compatible with the model. A local change of geometry or 'move' which 
is clearly ergodic (i.e. can generate any of the allowed configurations when applied 
successively) is shown in Fig. ^. It consists in deleting the two triangles adjacent 
to a given space- like link (if the resulting configuration is allowed). Its inverse is 
a 'split' of a given vertex and two neighbouring time-like links into two, thereby 
creating a new space-like link, as well as two new triangles. This is a special case 
of the so-called 'grand canonical move' sometimes used in dynamical triangulation 
simulations [112], O, 0], and does not preserve the total volume of space-time. 



Detailed balance equations for the move can be derived from standard consid- 
erations 0. Let Ny denote the number of vertices (A''^ = where is the 
number of triangles), and v a specific vertex. For pure Lorentzian gravity without 
matter, the equation for detailed balance reads 



P{Nv^Nv + l) 

Nykinkout 



P{Nv+l)P{Nv + l^Nv), 



(32) 



where P{Ny) = ^ is the probability distribution for labelled triangulations, 
and kin and kout count the incoming and outgoing time-like links at v (see Fig. 
1). We are still free to choose P{Ny ^ Ny^V) and P{Ny\ \ Ny) such that 
condition (|3^) is satisfied. Once a transition probability P(A'y + 1 — > A"y), say, has 
been chosen, it will be tested during the simulation against the uniform probability 
distribution between and 1 as follows. Choose a random number r g]0, 1]. Then, 
if the move is allowed (i.e. if the resulting triangulation belongs to the allowed class 
of configurations) it is accepted if P(A'y + l Ny) > r. If it is not allowed, one 
proceeds to the next move. 
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It is straightforward to generalize the updating of geometry to include Ising spins. 
The spin Hamiltonian is included in P{Nv), which now becomes a function of both 
Ny and the spin configurations. When inserting a vertex v, one has to specify at 
the same time a spin associated with v. The choice of spin up or down is made 
with probability 1/2, and the final result tested as in the case of the pure geometry 
update. 

We have performed the computer simulation for surfaces with toroidal topology 
and for system sizes of iV=2048, 4050, 8192, 16200 and 32768 triangles, and with 
a number t=32, 45, 64, 90 and 128 of time-slices respectively. Since the moves are 
not volume-preserving, fixing the system size to N is implemented as follows: we 
allow the volume to fluctuate within a certain, not too wide range, and collect for 
every sweep the first configuration with volume A^. The volume fluctuations are 
controlled by adding a term 6X (AN)"^ to the action, where AA^ is the deviation of 
the volume from its desired value A^. This term does not affect the ensemble of 
configurations collected, since for all of them AN = 0. We find that {AN)~^ ~ y/SX. 
Finally, one checks that the results obtained do not depend on the chosen, allowed 
range of volume fluctuations. A sweep is a set of approximately Ny accepted moves. 
For each /5-value used in the multi-histogramming analysis we perform 1.25 x 10^ 
sweeps (0.75 — 1.00 x 10^ for A^ = 32768). Measurements are made every 5 sweeps 
and errors are computed by data binning. 



3.1 Numerical results for the spin system 

The determination of the critical properties of the Ising spin system coupled to 
Lorentzian gravity proceeds in two steps (see for a recent, more complete dis- 
cussion in the context of 2d Euclidean quantum gravity). We first locate the critical 
/?-value where the system undergoes a transition from a magnetized (at large (3) 
to an unmagnetized phase. Next, we perform simulations in the neighbourhood of 
the critical value f5c and use finite-size scaling to determine the critical exponents. 
Finite-size scaling is also very useful for determining the location of the critical 
coupling itself, since a number of standard observables show a characteristic be- 
haviour for jS close to jSc- The following are some of the observables we have used, 
together with their expected finite-size behaviour (see for a full list): 



X = A^((m2) - {\m\f) ~ N^/'^'^" (susceptibility) (33) 

AnH=iv((e)-^)~iV""^» (A.H-^) (34) 

(....^), ,35) 

where 7 and u are the critical exponents of the susceptibility and of the divergent 
spin-spin correlation length, and dn is the Hausdorff or fractal dimension of space- 
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time. For a flat space-time (where of course dH = d = 2), we liave ud = 2, wliereas for 
the Ising model coupled to Euclidean quantum gravity udn = 3 (and dfj ^4). The 
internal energy density e and the magnetization m of the spin system are defined 

by 

NZn{P) dp ' NZn{P,H) dH H=Q ^ ' 

In order to find the critical point /5c, we can use the fact that the pseudo-critical 
coupling (5c{N) at volume N is expected to behave like 

/^^(^) - + (37) 

close to = Pc{N = oo), with c a constant. The observables (|33|)-(|35D all have 
well-defined peaks which we used for a precise location of /3c(iV), with the help of 
multi-histogram techniques. Eq. ( ^Tf ) can now be used to extract Pc and l/vdu- 
However, it is advantageous to determine first l/fdn from the peak values of (|34[) 
and (|35|) , and then substitute this value into (|37[) , thus reducing the number of free 
parameters in the fit. Afterwards, one can check that consistent values for l/vdu 
are obtained from the observables ( plD and ( pS]) at /3c, using multi-histogramming. 
In our simulations, l/vdu extracted from the peaks was so close to the Onsager 
value 1/2 that we did not hesitate to use l/z/(i/^ = l/2 in (pT)). Leaving it as a free 
parameter one obtains consistent results, but the error in j3c becomes larger. Finally, 
we have determined '^jvdu from (0), both from the peak values and at /3c (and 
again using multi-histogramming). 

In the table below we have listed the /3c- values extracted from measuring Z^inimi 
and Anm2, and assuming that Ijvdu = 1/2. For comparison with the high-T results, 
we also give the corresponding critical values for the expansion parameter u. Similar 
results are obtained from the rest of the observables we have measured. 



Observable 


/3c, using ^ = \ 


Uc = tanh /3c 


-Din |m| 


0.2522 (2) 


0.2470 (2) 


-Dlnm^ 


0.2520 (1) 


0.2468 (2) 



Column 1 of the following table contains the values of l/vdu extracted from the 
peak position for all three observables (|3^ ) -(|35| ) by using relation (|37|) (with free 
parameters /3c, l/vdu and c). Column 2 and 3 give l/vdn extracted directly from 
(P^) and (p5D by using the peak values of the observables and their values at (3c- 
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Observable 


^— from peak pos. 


-3— at peak 

vaH 


77777 (3c 


X 


0.52 (2) 






An|m| 


0.47 (2) 


0.531 (2) 


0.521 (3) 


^In 


0.53 (1) 


0.531 (1) 


0.520 (3) 



Lastly, we give the value of ^jvdu extracted from the susceptibility 



Observable 


value at peak 


value at (3c 


X 


0.883 (1) 


0.899 (2) 



Comparing the estimates for (3c from the high-T expansion and the Monte Carlo 
simulation, one finds good agreement. The results of the simulation strongly indicate 
that the critical exponents are given by the Onsager values vdu = 1/2 and 7 = 
1.75. Again, this corroborates the conclusion already reached by means of the high- 
T expansion. Further evidence that the system belongs to the Onsager and not 
the Euclidean gravity universality class comes from measuring the magnetization 
exponent (3m/i'dH and the specific heat exponent ajvdu- Their Onsager values are 
1/16 and 0, whereas in Euclidean gravity they are 1/6 and —1/3. In our model, 
the magnetization exponent determined from ~ j^-Pmlydn found to 

be (3mlydu = 0.070(1), favouring the Onsager value 0.0625 over the Euclidean 
gravity value 0.1666. The specific heat exponent was obtained from the finite-size 
scaling of the values of the specific heat peaks Cy ~ ]\joi/udH _ power fit yields 
a/vdu = 0.0861(7) at x^/dof = 11.6 whereas a logarithmic fit gives x^/dof = 1.57, 
supporting the conjecture that a = 0. 

We do not have independent measurements of the critical parameters v and dn 
from the spin sector alone, but we will determine the Hausdorff dimension du in 
the next subsection from an analysis of the geometry of Lorentzian quantum gravity 
coupled to Ising spins. 



3.2 Numerical results for the geometry 

As is well-known from both analytical studies |15|, O, |19[ and numerous Monte Carlo 



simulations (^171 16 



and references in g]), finite-size scaling is a powerful tool for 
determining the fractal space-time structure of two-dimensional Euclidean quantum 
gravity. The same technique can be used to investigate the geometric properties of 
two-dimensional Lorentzian quantum gravity. 

In a given triangulation, we define the distance between two vertices as the 
minimal length of a connected path of links between them. In 2d Euclidean quantum 
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gravity this notion of distance becomes proportional to the true geodesic distance 
between the vertices in the infinite- volume limit. We will assume this is also true 
for the present model. All diffeomorphism-invariant correlation functions of matter 
fields must be functions of this geodesic distance. Both the geodesic distance and 
the fractal dimension appear in the expression for the volume 

Ar(r)~r'^« for r < A^^/'^«, (38) 

which denotes the number of vertices (or triangles) inside a ball (a disc in dimension 
2) of link-radius r. If n^{r) denotes the number of vertices at distance r from a fixed 
vertex v, relation ( |5BD implies that 

ny{r) r'^"-^ for r < iV^/'^«. (39) 

Finite-size scaling for an observable A then leads to a scaling of the correlation 
function integrated over all points at distance r from a vertex according to 

{A{r)A{0))j^ N^-'/'^H-^^FAix), x = r/N^/'^». (40) 

The factor N^~^I'^h comes from the integration over points at distance r from vertex 
using (|39|), while A^^ is the genuine dynamical exponent of the correlator. 
By measuring correlation functions for various volumes iV, one can determine 
dfj and the critical exponents. We will concentrate here on the Hausdorff dimension 
dij. One first rescales the height of the measured distributions {A{r)A{Q))^ to a 
common value. However, the distributions measured for different N will still have 
different width as functions of r. By appropriately rescaling r, they can then be 
made to overlap in a single, universal function Fa{x). From a technical point of 
view it is important to work with the shifted variable 

X = (41) 



where the shift may depend on the observable A. Using eq. (^Tj) takes into account 
in an efficient way the major part of the short-distance lattice artifacts, as has been 
discussed carefully in [jl8|, |19| . Applying standard procedures from Euclidean 2d 



quantum gravity then leads to the results summarized in Table ^ The observables 
appearing in Table ^are: (i) the number n„(r) of vertices at a given (link-) distance 
r from a fixed vertex v, which may be viewed as the correlation function of the 
unit operator in quantum gravity [|l^]; (ii) the number Sup{r) of spins at distance 
r from a vertex v which are aligned with the spin at f ; (iii) the number Srfo^„(r) 
of spins with orientation opposite to the spin at v; (iv) the spin-spin correlation 
function s{r) between vertices separated by a geodesic distance r; (v) the function 
S{r), obtained by integrating s(r) over all vertices at distance r from a vertex 
(vi) the distribution SV{1) of spatial volumes, with / denoting the length of a given 
time-slice. For the shift a^, we obtained the estimate — 4 < < — 1. Unfortunately 
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Figure 6: The distribution of spatial volumes SV{1) at [3 = [3c, rescaled according to 
eq. (13). We have set dH = 2.0. 



our statistics was not good enough to determine it with more accuracy. However, 
the fact that it is non-vanishing justifies its introduction in the first place. After a 
suitable normalization, we expect the volume distribution to behave like 

SV{1) ~ fil/N^''^"). (42) 

Fig. I demonstrates clearly that for the Ising model at /? = /?c, SV{1) scales as 
anticipated when we set (i// = 2.0. Scaling the Ising distributions at [3 = (3c{N), the 
pseudo-critical point of the magnetic susceptibility, or considering pure gravity leads 
to similar results. 

We conclude from Table |l| that the Hausdorff dimension of 2d Lorentzian quan- 
tum gravity is close to the flat-space value dH = 2. This is clearly different from the 
Euclidean situation which is characterized by = 4 for pure gravity and c?// > 4 in 
the presence of Ising spins. The results for the Lorentzian gravity-matter system are 
particularly convincing for the purely geometric observables n^{r) and SV{1), which 
basically coincide with the corresponding measurements obtained in Lorentzian pure 
gravity. 



4 Conclusions 

We have presented compelling evidence, coming from a high-T expansion as well 
as Monte Carlo simulations, that the critical exponents of the Ising model cou- 
pled to Lorentzian gravity are identical to the exponents in flat space. This is in 
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Observable 


C 2 ' ''^c 


dn 

c = i, /3 = ^{N) 


c = 




2.00(5) 


2.00(5) 


2.03(4) 




1.92(5) 


2.08(2) 


— 


^down 


2.20(3) 


2.07(2) 




s 


2.10(3) 


2.20(5) 




s 


2.05(7) 


2.05(5) 




sv 


2.00(4) 


2.00(3) 


2.00(3) 



Table 1: The Hausdorff dimension dn, obtained from the Ising model scaling at (3^ 
(column 1), at (3c{N) (column 2) and from pure gravity (column 3). 

contrast with the situation in Euclidean gravity (i.e. Liouville gravity), where the 
exponents change^. Similarly, the fractal dimension of space-time is unchanged in 
the Lorentzian model after coupling it to a conformal field theory (the Ising model 
at the critical point). In Euclidean gravity the fractal properties of space-time are 
in general a function of the central change of the conformal field theory. From the 
evidence collected so far, we conclude that matter and geometry couple weakly in 
Lorentzian gravity and strongly in Euclidean gravity. 

For the case of the Ising model, this difference can be explained in more detail in 
geometric terms. As mentioned earlier, it has been shown in [|I| that the difference 
between Euclidean and Lorentzian gravity is related to the presence or absence 
of baby universes. On the other hand, it is by now well understood that baby 
universes are at the source of the strong coupling between spins and geometry. This 
can happen because the spin configuration of a baby universe can be flipped relative 
to that of the 'parent' universe at almost no cost in energy since the 'baby' and 
the 'parent' are connected only by a thin tube. The geometry of two-dimensional 
Euclidean quantum gravity is very fractal, with many 'pinches' at all scales, leading 
to typical spin configurations that look very different from those on fiat space-time. 
Moreover, the presence of Ising spins on the surface effectively enhances the fractal 
baby universe structure since it is exactly the lowest energy spin configurations 
(apart from the ground state) that involve baby universes. The interaction becomes 

^The exponents of the Ising model coupled to 2d Euclidean quantum gravity are equal to those 
of the 3d spherical model. It is not understood whether this is a coincidence. More generally, the 
exponents of the (m, m + 1) minimal conformal model coupled to 2d Euclidean quantum gravity 
agree with the critical exponents of the spherical model in 2m/ (m +1) dimensions. 
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so strong that it tears the surface apart when we couple more than two Ising spins 
to the two-dimensional geometry. This is the origin of the famous c = 1 barrier of 
two-dimensional Euclidean quantum gravity. 

Once the creation of baby universes is disallowed, as in the case of the Lorentzian 
model, the coupling between matter and geometry becomes weak, and the matter 
theory has the same critical exponents as in flat space-time. This happens although 
the typical space-time geometry is by no means fiat, a fact we have already empha- 
sized in the introduction, and which is illustrated by Fig. |I|. On the contrary, our 
model allows for maximal fluctuations of the spatial volume which can jump from 
(essentially) zero to infinity in a single time step. However, as we have been able 
to demonstrate, such violent fluctuations of the two-geometry are still not sufficient 
to induce a change in the critical exponents of the Ising model. From the above 
arguments it seems likely that Lorentzian gravity can avoid the c = 1 barrier. This 
question is presently under investigation. 
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